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ABSTRACT 

We cross-correlate the cosmic microwave background temperature anisotropy maps from the WMAP, 
MAXIMA-1, and MAXIMA-II experiments. We use the cross-spectrum, which is the spherical har- 
monic transform of the angular two-point correlation function, to quantify the correlation as a function 
of angular scale. We find that the three possible pairs of cross-spectra are in close agreement with 
each other and with the power spectra of the individual maps. The probability that there is no 
correlation between the maps is smaller than 1 x 10^^. We also calculate power spectra for maps 
made of differences between pairs of maps, and show that they are consistent with no signal. The 
results conclusively show that the three experiments not only display the same statistical properties 
of the CMB anisotropy, but also detect the same features wherever the observed sky areas overlap. 
We conclude that the contribution of systematic errors to these maps is negligible and that MAXIMA 
and WMAP have accurately mapped the cosmic microwave background anisotropy. 
Subject headings: cosmology: cosmic microwave background, methods: statistical, methods: data 
analysis 



1. INTRODUCTION 

Temperature fluctuations in the cosmic microwave 
background (CMB) encode a vast amount of cosmo- 
logical information about our universe. CMB pho- 
tons released from the primordial plasma at the time 
of recombination approximately 380,000 years after the 
Big Bang, provide thus a picture of the universe in 
its infancy only somewhat modified by low-redshift ef- 
fects such as reionization. Recently WMAP produced 
a 13' full sky measurement of the CMB tempera- 
ture anisotropy ({Bennett et alJl2003a() . This map has 
been used in conjunction with other CMB and cos- 
mological data to constrain a number of cosmologi- 
cal parameters to unprecedented accuracy ijSpergel et al.l 
I2OO3') . Previous to WMAP a number of exper- 
iments produced high quality maps of CMB tem- 
perature anisotropy. These included both balloon 
borne bolometric experiments such as B OOMERANG 
(|de Bernardis et al."'2000'; ''Ruhl et a l (m)2l) , MAXIMA- 
I fflanany et a l. 2000; L ee et al. ,200lj) . and ARCHEOPS 
(jBenoit et all l2003l). and ground based interferqmetr ic 
experiments, CBI (iPadin et alJl2001l:rMason et alJl200311 
DASI i Halvcrs oiTet alJ I2002D and VSA ijGrainge et all 
l2002l|) . Tight constraints were placed on cosmological pa- 
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rameters from these experiments as well (e.g. [J affc ct alJ 
Imil: iNetter field et al., (,.2002)_: Stomnor et al. (200L); 
lAbroe et alJ l|2002D: IPrvke et aTl (ppavT 

Given the longer observations and higher sensitivities 
of recent CMB experiments the theoretical analysis is 
more likely to be limited by systematic rather than sta- 
tistical errors. It is therefore important to ensure that 
systematic errors in these experiments are sub-dominant 
compared to statistical errors. A comparison of the 
power spectra from the experiments can provide some 
confidence that systematic errors are not dominant. For 
instance the power spectra of WMAP, MAXIMA-I, and 
MAXIMA-II are shown in Figure [H Note that no cal- 
ibration adjustments have been made to the data. For 
experiments which observe overlapping parts of the sky 
the close agreement of the power spectra does not nec- 
essarily imply that the spatial fluctuations detected by 
the experiments are identical. In such cases, and partic- 
ularly when the experiments have similar angular resolu- 
tion, a more stringent test for systematic errors is to cross 
correlate the temperature fluctuations of one map with 
the fluctuations in the other. Positive correlations be- 
tween temperature anisotropy maps would also enhance 
the confldence in the reconstruction of the spatial pat- 
tern of the CMB. Some difficulties arise if the two ex- 
periments under consideration have different pixel reso- 
lutions and beam profllcs, which is usually the case. In 
that case a straightforward pixel to pixel comparison is 
no longer accurate because the CMB signal contained in 
corresponding pixels is not the same, and a more elabo- 
rate technique needs to be employed. 

In this paper we use the cross-spectrum as a technique 
to cross-correlate the maps of WMAP and MAXIMA. 
This is the flrst reported cross-correlation of a CMB map 
with the map of WMAP, and the flrst release of data 
from the MAXIMA- JJ flight. 

Other comparisons of CMB anisotropy maps 
from different experiments have been performed in 
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Fig. 1.— The CMB angular power spectra as measured by WMAP fHinshaw et al.' ^2003|), black circles), MAXIMA-I ILee et alj 120011) . 
blue squares), and MAXIMA-II (red diamonds). No adjustments have been made to the calibration for any of the spectra. For the analysis 
in this paper we use the version of the MAXIMA -1 data published bv nHananv et al.. 1.20001) (see Section y}- 



the pa st (e.g. iGanga et alJ l)1993() . Tenerife and 
DMR llLinweaver et alJ \199m . MSAM and Saska- 
toon fKnox et al.' '199g), and QMAP and Saskatoon 
ijXu. T cgmark. & dc 01iveira:£psta 200^). However, 
none of these analyses used the cross spectrum as 
a technique for quantifying the amount correlation 
between data sets. 

This paper is organized as follows: in Section[21we give 
details of the cross spectrum and its use for quantifying 
the amount of correlations between two CMB anisotropy 
maps at various angular scales. In Section Owe discuss 
the maps used in our analysis. Results of the calcula- 
tion of the power and cross-spectra for the maps, and 
an analysis of the difference maps are given in Section^ 
We discuss the results in Section |S1 and summarize in 
Sectional 

2. METHOD 

Perhaps the simplest way to compare quantitatively 
two maps is to calculate a for their difference. This 
statistic has the limiting feature that only overlapping 
sections of the maps can be included, and the maps must 
have identical pixelizations. Additionally, if the maps 
have different beam window functions, then the signal in 
corresponding pixels is different and the distribution of 
the statistic would no longer be distributed. The 
statistical interpretation of the value would therefore 
be difficult. 

There are a number of statistics which can be calcu- 
lated between two CMB maps to quantify their consis- 
tency, e.g. the li near and rank correlation coefficients 
ijPress et al.lll99^ . Unfortunately, these correlation co- 
efficients also suffer from the same difficulties as the 
statistic. 

Several authors derived statistics which do take into 



account partial overlapping maps, diff e rent p ixeliza- 
tions and beam profiles. iKnox et al.l 1)19981) derive 
both Bayesian and frequentist techniques for correlating 
CMB maps, and apply these statistics to data from the 
MSAM92, MSAM94, and Saskatoon experiments. They 
advocate the calculation of the "contamination parame- 
ter" , which gives the probability distribution for the mag- 
nitude of a signal that is not common between the two 
data sets under consideration. A low value for the con- 
taminat ion parameter im plies that the data sets are con- 
sistent. iTegmarfl l)1999|l defines a "null-buster" statistic 
which gives the number of a between the difference of 
two maps and a hypothesis of pure noise. These statis- 
tics can be used both as an internal consis tency check be- 
tween detectors for the same experiment (jStompor et alJ 
and to compare maps from different experiments. 

In this paper we use the cross-spectrum to quantify 
the level of correlations between the data of WMAP and 
MAXIMA . The cross-spectrum is the spherical harmonic 
transform of the real space correlation function for the 
two maps. Whereas other statistics condense the in- 
formation about correlations into a single number, and 
therefore result in some loss of information, the cross- 
spectrum retains more information by analyzing the cor- 
relations as a function of angular scale. 

We now discuss the method for estimating the cross- 
spectrum from CMB temperature anisotropy maps. 
Consider two maps of the CMB called T^^' and T'^), 
respectively. Then 

AT;W=,sW+n|^) (1) 
Ai;(^)=,sf)+nr\ (2) 

where j is a pixel index, and s.p'' are the CMB signal 
in pixel i, and n^''"'' and are the pixel noise for the 
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first and second map, respectively. Let the number of 

pixels in the first and second maps be A^p''"'' and nI^\ 
respectively. We write the data vector in pixel space as 



(3) 



(1) 



(2) 



where d is now a column vector of length Np 
Assuming the signal and noise within each experiment 
are uncorrelated and that the noise between experiments 
is uncorrelated, and using Equationsnand|21we find that 



{d(f)^M = 



g{C)T 



5(2)+_^(2) 



(4) 



where AA(i) = {n^^^n^^^^) andAA(2) = (n(2),i(2)T^ 
pixel noise covariance matrices for the first and second 
experiment, respectively. The quantities S^^\ S^'^\ and 
S'(*-') are the CMB signal covariance matrices, which can 
be written as 
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-Cf^i3l')i?f'p,(cos%), (7) 



auto-spectrum estimated for any one experiment alone 
may differ from the auto-spectrum calculated when es- 
timating the cross spectrum. For the case of a full sky 
coverage one can show that the level of correlations be- 
tween the spectra is proportional to the amplitude of the 
cross-s pectrum, that is to the am ount of common sky 
signal l)Kamionkowski et alJll996() . We therefore expect 
a high level of correlations between the spectra if they 
share the same sky signal at low fs where the contribu- 
tion of instrument noise is smaller compared to sample 
variance. 



MAXIMA- 1 & 
MAXIMA -n 



where b''^^ and Bf^ are the beam profiles for the first 
and second experiment, are Legendre polynomials, 
and 6ij is the angle between pixels i and j. The auto- 



spectra c'}^'^ and C^"' are commonly called the power 
spectra for the first and second map, respectively, and 
C\ is defined as the cross-spectrum. In the zero noise 
case the cross-spectrum is limited by the requirement 

Otherwise M would have neg- 
ative eigenvalues and therefore be unphysical. Up to an 
irrelevant additive constant, the likelihood for d is 

- 2 \nC={N^^^ + iV^2)) In \M\ + d^Ar^d. (8) 

We maximize the likelihood £ as a function of the 
two auto-spectra and the cross-spectrum. If two CMB 
maps have a high degree of correlation, then we expect 
the cross-spectrum to resemble the auto-spectra of both 
maps. If the signal in the two maps is not correlated 
the cross-spectrum should be consistent with zero at all 
angular scales. 

To maximize C we adopt a Newton-Raphson technique 
for finding the zero of the first derivative of the log like- 
lihood funct ion. We use a technique sim ilar to the one 
described bv lBond. Jaffe. fc Knoxl l)1998() except that we 
use the full curvature matrix instead of the Fisher ma- 
trix when calculating the s teps for convergence of the 
Newton-Raphson algorithm (|Hobson fc M aisingcr 2002). 
We found that use of the Fisher matrix gave non-positive 
definite pixel-correlation, presumably because the likeli- 
hood function of the three spectra together has a com- 
plicated structure. 

The cross-spectrum method involves estimating all 
three power spectra simultaneously and gives rise to cor- 
relations between the different spectra. Therefore an 
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Fig. 2.— The beam filter functions for MAXIM A-I (dotted line), 
MAXIMA-II (dashed line), and WMAP (solid line) considered in 
this analysis. The horizontal axis is spherical harmonic multipole 
as a function of £. 



The cross-spectrum is a powerful technique for com- 
paring two CMB anisotropy maps because it accounts 
for different beam shapes, pixel resolutions and sky cov- 
erage in a simple and straightforward way. For exam- 
ple, one can compute the cross-spectrum for two maps 
that do not overlap at all. In such a case the results 
would not be sensitive to correlations on small angular 
scales. Also, though our formalism describes estimating 
the cross-spectrum for two maps, a further generalization 
could be made to estimating correlations between three 
or more maps simultaneously. However, given numeri- 
cal subtleties which are subsequently discussed in Sec- 
tion we found it prudent, while equally convincing, 
to estimate the correlations for only pairs of maps. 

The power spectra estimated from the WMAP maps 
are actually the cross-spectra fr om different de tectors 
from the same frequency band ijHinshaw et al.l |2003) . 
Our approach is distinct from the one used by the 
WMAP team in t hat they use a frequentist approach 
(jHivon et alJ200^ to estimate the cross-spectra between 
various detectors of the same frequency band indepen- 
dently from the auto-spectra. The approach we use in 
this paper is entirely Bayesian. 

3. THE MAPS 

We cross-correlate CMB temperature anisotropy maps 
from the WMAP, MAXIMA-I and MAXIMA-II experi- 
ments. The WMAP map used is the W-band (93 GHz) 
foreground-cleaned map, and only the portion that over- 
laps with the MAXIMA-I field. By analyzing 





Fig. 3. — A comparison of the overlap region of the MAXIMA-I (upper right panel) and the WMAP 93 GHz (upper left panel) CMB 
temperature anisotropy maps, and their difference map (lower panel). The maps have been Wiener filtered for a visual comparison; the 
non-filtered versions are used in the analysis. All modes with ^ < 35 have been removed from these maps. The MAXIMA-I map, which 
originally had a resolution of 10', has been smoothed to the WMAP resolution. The color scale is from —250 to 250 fiK. 



WMAP 's derived foreground maps we find that this por- 
tion of the W-band map is free of point sources, and 
contains a neghgible amount of dust, free-free, and syn- 
chrotron emission^. The WMAP data at 93 GHz has 
similar angular resolution and is closest in frequency to 
the MAXIM A data. This map is pixelized using the 
HEALPixio l|Gorski. Hivon. fc Wandell IT999(l pixehza- 
tion in Celestial coordinates with nside= 512, and con- 
tains 7,926 pixels. An nside=512 corresponds r oughly 
to a pi xel with a width of 7'. Following Benne tt et al.l 
^^2003a^ we assume that there are no noise correla- 
tions between pixels and that the beam pattern is non- 
Gaussian, azimuthally symmetric and has FWHM of 
13.2'. We computed the beam window function by tak- 
ing the weighted average of the individual beams for each 

^ http:/ /lambda. gsfc.nasa.gov/product/map/m_products.cfm 
http:/ /www. eso.org/science/healpix/ 



of the four detectors used to form the final W-band map. 
The resulting beam profile in £ space is shown in Fig- 
ure 121 The Wiener filtered version of this map is shown 
in the left panel in Figure |3 we use the raw data in the 
analysis. The Wiener filtering is performed with the cor- 
responding best fit models for each map. A more detailed 
discuss ion of the WMAP maps is given in lBennett et alJ 
l2003aD and of algori thms used in their computations in 
iHinshaw et"an l)2003|) . 

The MAXI MA map-making proce dure is described ex- 
haustively in iStompor et"all l|2002D . The MAXIMA-I 
map we use in this analys is is the 8 arc r ninute version 
of the data published by iHananv et al.l 1)2000(1 . which 
covers a larger area of the sky and has a coars er res- 
olutio n compared to the data published by iLee et alJ 
ll200l . There are a total of 5,972 pixels, and this map 
covers ~ 100 deg^ on the sky. It is a combination of 
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Fig. 4. — The overlap sections of the MAXIMA-I and MAXIMA-II Wiener filtered maps, shown in the left and center panel, respectively. 
The difference of the Wiener filtered maps is shown in the right panel. The color scales range from — 300/.tK to 300^K. Note that while 
only the overlap section of the MAXIMA-I map is shown here, we use the entire MAXIMA-I map for all analyses. 



three 150 GHz photometers and one 240 GHz photome- 
ter. The beams for ea ch photometer have a FWHM 
of ~ 10' ijHananv et al .1 ^2000') and the cfTective window 
function i s calcu lated using the technique described in 
iWu et alJ 1)200 H) . The MAXIMA-I map shown in Fig- 
ure is Wiener filtered and smoothed to a WMAP reso- 
lution. The raw version of the map is used for all quan- 
titative analyses. Figure |3| also shows the difference of 
the MAXIMA-I and WMAP Wiener filtered maps. The 
pattern of temperature fiuctuations, which is similar in 
both maps, disappears in the difference map. 

The MAXIMA-II map comes from a flight of the MAX- 
IMA payload that took place on 1999 June 17. We use 
the data from four photometers at 150 GHz, and only 
the portion of the MAXIMA-II map that overlaps the 
map of MAXIMA-I More details about the MAXIMA- 
II flight, data and maps are given in Rabii et al. (2003, in 
preparation) , and'S tompor et all H2003l|l . The MAXIMA- 
II map is pixelized using an 8' square pixelization in ce- 
lestial coordinates, contains 2,757 pixels, and covers ~ 
50 deg^ on the sky. The beam profile for this map is 
~ 10' FWH M, and aga i n com puted using the techniques 
described in 'Wu et al.' (2001). The MAXIMA-II power 
spectrum shown in Figure ^ has 10 bins of = 75, ex- 
tending over the I range 35 < £ < 785. Figure^ shows 
the overlap region of the MAXIMA-I and MAXIMA-II 
maps and the difference map. Identical temperature fiuc- 
tuations that are apparent in each of the maps disappear 
in the difference map. 

4. RESULTS 

The auto- and cross-spectra for all combinations of the 
WMAP and MAXIMA maps are shown in Figure 13 The 
error bars on the spectra are the square root of the cur- 
vature of the likelihood function about the maximum 
likelihood parameter value. In all cases we compute 
the spectra in bins of width A£ = 75, over the interval 



Table 1 
Cross Spectrum Values 



Maps 


DoF 


Xcutoff 




MAXIMA-I /WMAP 


8 


53 


191 


MAXIMA-I /MAXIMA-II 


8 


53 


241 


MAXIMA-II /WMAP 


8 


53 


150 



The from EquationElcalculated for all three combinations 
of maps. A greater than 53 implies that the probability 
that the no correlation hypothesis is true is less than 1 x 10~*. 



Ill < ^ < 710, and marginalize over all modes ^ < 110 
and (i> 711. The appropriate pixel window functions for 
each map were convolved with the beam functions in the 
analysis. We found that the cross-spectrum estimator 
did not converge when the initial bin was split in two, 
and this is further discussed in Section|Sl In all cases the 
cross-spectra are consistent with the auto-spectra giving 
strong evidence for a correlation between the maps. 

We also compute the power spectrum of the difference 
maps for all three pairs of maps using bins of A£ — 75 
over the range 35 < ^ < 785. The WMAP window func- 
tion is used when computing the MAXIMA-I - WMAP 
and MAXIMA-II - WMAP difference spectra. The ex- 
pected residual power resulting from different beam pro- 
files is maximum at the bin centered at £ ~ 300, and is 
approximately equal to the la error bar of the MAXIMA- 
I / WMAP difference power spectrum. The effect is less 
than 1(7 for all remaining bins. We use the MAXIMA- 
I window function when computing the MAXIMA-I - 
MAXIMA-II difference spectrum. The results are shown 
in Figure Of the 30 band power estimates for the 
difference maps, 28 are within Icr of zero power. 
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Table 2 

Difference Spectrum Values 



Maps 


DoF 




MAXIMA-I /WMAP 


10 


7.5 


MAXIMA-I /MAXIMA-II 


10 


8.3 


MAXIM A-II /WMAP 


10 


17.2 



The of the power spectrum for the difference maps with 
the null spectrum. 



To further quantify the level of correlation between the 
maps we use a statistic to reject the hypothesis that 
the maps are uncorrelated. We write our statistic as 

X^ = E^b"^^--^b"^ (9) 

BB' 

where the sum is over band power estimates, and F is 
the Fisher matrix for the cross-spectrum. Because the 
auto-spectra and cross-spectrum are estimated simulta- 
neously, we marginalize over the auto spectra when cal- 
culating the x^- 

To test the null hypothesis we choose a statistical sig- 
nificance a = 1 X 10~*. If is greater than the critical 
value 53, then the probability that the null hypothesis is 
true is less than 1 x 10~*. The results are summarized 
in Table ^ In all cases is significantly larger than 
the critical value giving an essential certainty that the 
no-correlation hypothesis is false. Note that assuming 
the x^ of Equation |H| is x^ distributed is equivalent to 

(C) 

assuming that the C^' are Gaussian distributed, which 
is an approximation. 

We also compute the x^ of the difference spectra shown 
in Figure with the null spectra to determine how con- 
sistent these are with no fluctuations in the difference 
maps. The results are shown in Table 12 The 10 power 
spectrum bins computed from the MAXIMA-I /WMAP 
and MAXIMA-I /MAXIMA-II difference maps have a 
of 7.5 and 8.3, respectively, with the null spectrum. The 
MAXIMA-II /WMAP difference map gives a x^ of 17.2. 
There is a 7% chance of getting x^ > 17.2 for 10 DoF. 
Overall there is a good fit to the null spectrum model, 
which implies that differencing the overlap section of the 
maps removes the sky signal and is consistent with noise. 

5. DISCUSSION 
5.1. Auto- and Cross-Spectra 

The auto- and cross-spectra of the different data sets 
agree with each other to within la over almost all £ bins 
giving evidence that at each angular scale all experiments 
are detecting the same spatial fluctuations on the sky. 
All auto- and cross-spectra show the first acoustic peak 
in the power spectrum and then a level of power that 
is consistent with subsequent peaks. These results are 
consistent with standard inflationary ACDM models. 

Auto-spectra of the overlap section of the WMAP data 
give increased error bars at £ > 486 because of the limited 
sky coverage of the overlap regions, and because of the 
beam proflle of the W-band map. We flnd that the beam 
pattern alone causes the WMAP auto-spectrum error 
bars in the bins ^={{486, 560}, {561, 635}, {636, 710}} to 



be 2-3 larger than those for MAXIMA-I or MAXIMA-II. 
Negative power was found in the bin ^={486, 560} for 
the WMAP auto-spectra (see the top and bottom panel 
of There is no requirement that the auto-spectrum 
be positive in our estimation method. 

A comparison of the auto-spectra shown in Figure El 
reveals that there is a difference between band power 
estimates for the same dataset. This difference arises 
because the computation of cross-spectra involves es- 
timating both auto- and cross-spectra simultaneously, 
giving rise to correlations between the different spec- 
tra. The fractional changes in power averaged over 
bins are 3%, 6%, and 15% for the WMAP, MAXIMA-I, 
and MAXIMA-II data sets, respectively. If the likeli- 
hood distribution of the band powers were strictly Gaus- 
sian, then the maximum likelihood estimates would be 
the same regardless of the correlations between spec- 
tra. However, the likelihood as a function of auto- and 
cross-spectra is somewhat n on-Gaussian (see Equation|Hl 
IBond. Jaffe. fc Knoxl ^MB)), so the correlations do ef- 
fect the band power estimates. The fact that the changes 
between estimates are small suggests, however, that the 
distributions are close to Gaussian. 

We note that the bin-powers are correlated at the 77% 
level or higher at the lowest i bin. The correlation de- 
creases to less than 10% at the highest £ bin for corre- 
lations between the auto-spectra and to 20% - 50% for 
correlations between the auto- and cross-spectra. These 
results are in broad agreement with expectations given 
the high amplitude of the cross-spectrum. 

5.2. Computational Issues 

The strong correlation between the different data sets 
leads to some computational difficulties when attempting 
to find the maximum likelihood auto- and cross-power 
spectra. As discussed in Section |21 the cross-spectrum 

is limited by the requirement that C^'^^' < \J C^^^C^^\ 
We find that using a quadratic estimator (Fisher matrix) 
method for calculating the Newton-Raphson step leads 
to a mis-estimate of the step for SC^ when starting 
with a guess significantly far away from the peak in like- 
lihood space. This consequently results in a step which 
leads to a non-positive definite pixel covariance matrix. 
For example, an initial guess of a null spectrum leads to 
an unphysical pixel covariance matrix in all three cases 
we are considering. This is remedied by using the curva- 

ture matrix to compute SC^ , as discussed in Section |21 
Once the parameter values become sufficiently close to 
the maximum likelihood values either the curvature ma- 
trix or flsher matrix can be used to flnd the maximum 
likelihood parameters. Both techniques converge to the 
same set of parameters for all three analyses. 

The power spectra shown in Figure have been cal- 
culated with a broad initial bin at ^ = {2,110}. This 
is because the Newton-Raphson likelihood maximization 
technique did not converge if this low £ bin was split 
to two. Some binning structures would cause negative 
eigenvalues in the curvature matrix of the parameters 
or steps in parameter space that would lead to a non- 
positive dcflnite pixel covariance matrix, both of which 
are unphysical. 

We carried out simulations and found a similar phe- 
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Fig. 5. — The auto-spectra and cross-spectrum estimated for the 
MAXIMA and WMAP maps. From top to bottom: the results for 
MAXIMA-I /WMAP, MAXIMA-I /MAXIMA-II, and MAXIMA- 
II /WMAP . Note that only the portion of the WMAP map which 
overlaps with the MAXIMA-I field is used in the analysis. 



nomenon. The cross-spectra of uncorrelated maps or 
maps with a smaU value for the ratio of expected cross- 
spectrum to auto-spectrum converged to the expected 
answer. The calculation of the cross-spectrum also con- 
verged with simulated maps that had perfect correlation 
(i.e. the same map with different noise realizations) and 



Fig. 6. — The power spectrum of the difference maps from 
all three data sets. The diamonds (red) are the power spectrum of 
the MAXIMA-I / WMAP difference map, the squares (blue) are the 
power spectrum of the MAXIMA-I /MAXIMA-II difference map, 
and the triangles (green) are the power spectrum of the MAXIMA- 
II /WMAP difference map. All three spectra are statistically con- 
sistent with the null spectrum. 



a broad first bin with £ — {2, 110}. However it did not 
converge with simulated maps that had perfect correla- 
tion and two bins between ^ of 2 and 110. Therefore, 
we attribute the computation problems encountered as a 
limitation in the method used for computing the cross- 
spectrum and not a feature in any of the data sets con- 
sidered in this analysis. 

5.3. Foregrounds and Systematic Errors 

The cross spectrum and difference spectrum analyses 
conclusively demonstrate that all three experiments have 
mapped the same temperature fluctuations on the sky. 
However, these analyses are not sensitive to whether the 
shared fluctuations are CMB in origin or the result of 
foreground contamination or a shared systematic effect. 

A careful foreground analysis was carried out by both 
the WMAP and MAXIMA teams. It was shown that 
the MAXIMA-I region of the sky at 150-240 GHz con- 
tains a negligible amount galactic contaminatio n and 
that it has no detecta ble point sources IHanan v et alJ 

EHoS II^!iIir"2oo3) . 

A detailed analysis of th e foreground sources i n the 
WMAP data is presented in lBennett et alJ l)2003bD . Al- 
though we use WMAP 's foreground-cleaned map, 
even WMAP 's foreground maps have negligible amount 
of contamination in the MAXIMA-I region. The RMS 
fluctuations in the WMAP 93 GHz dust, synchrotron, 
and free-free maps are lower than the corresponding 
WMAP CMB map by a factor of 17, 42, and 560, re- 
spectively. Also, the WMAP team find no point sources 
in the MAXIMA-I region of the sky. 

It is unlikely that WMAP and MAXIMA share sys- 
tematic errors. We therefore conclude that the common 
signal in the WMAP and MAXIMA data is the cosmic 
microwave background radiation. Since the cross-spectra 
agree with the auto-spectra we conclude that within the 
signal-to-noise ratio of the tests systematic errors in the 
data are smaller compared to statistical errors. 
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6. SUMMARY 

We have presented a Bayesian method for estimat- 
ing the cross-spectrum between two CMB temperature 
anisotropy maps. The method is advantageous for corre- 
lating maps because it does not require the maps to have 
perfect overlap, identical beam shapes or pixelizations. 
Using this formalism we found a high degree of correla- 
tion between the maps from MAXIMA-I, WMAP, and 
MAXIMA-II ; in all cases the null hypothesis is rejected 
with a probability higher than 1 — 10~^. Additionally, we 
computed the power spectrum of the difference maps for 
all combinations of the three data sets considered, and 
found that in each case the spectra were consistent with 
the null spectrum. 

The results show conclusively that the tempera- 
ture fluctuations detected by each of the MAXIMA-I, 
WMAP, and MAXIMA-II experiments are reproduced 
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